Spin-orbit coupling and spin transport 
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Recent achievements in semiconductor spintronics are discussed. Special attention is paid to spin- 
orbit interaction, coupling of electron spins to external electric fields, and spin transport in media 
with spin-orbit coupling, including the mechanisms of spin-Hall effect. Importance of spin-transport 
parameters at spin-precession wave vector k so is emphasized, and existence of an universal relation 
between spin currents and spin accumulation at the spatial scale of l so w k^ 1 is conjectured. 
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INTRODUCTION 

Spin is the only internal degree of freedom of electron. 
Employing it for creating electrical and optoelectronic 
devices with new functionalities is the ambitious goal of 
spintronics . It is non less important that novel physics 
including concepts that are exciting and challenging but 
sometimes controversial are involved. Last several years 
have witnessed remarkable experimental and theoretical 
achievements in this field. 

Operating spintronic devices includes injecting 
nonequilibrium spins and manipulating spin polarization 
at given locations. Many of original proposals included 
electrical spin injection from magnetic (metallic or 
semiconductor) electrodes, similarly to spin injection 
into paramagnetic metals [2j, and manipulating spin po- 
larization by ac magnetic fields via Zeeman interaction. 
Modern trends in semiconductor spintronics are based 
on employing spin-orbit (SO) coupling for achieving 
both goals. It is believed, that such devices will provide 
efficient spin injection, and will be (i) controllable 
by electric fields allowing access to electron spins at 
nanometer scale, (ii) free from the parasitic effect of 
stray magnetic fields, and (iii) compatible with existing 
semiconductor technologies. 

With the SO coupling as the central paradigm of semi- 
conductor spintronics, developing a consistent theory of 
spin transport in systems with SO interaction is a de- 
manding goal. The history of a related phenomenon, 
the anomalous Hall effect (AHE), indicates existence of 
fundamental theoretical problems. As different from the 
regular Hall effect, the AHE is driven by electron mag- 
netization M rather than by an external magnetic field 
B , and therefore is based on SO coupling. The theory of 
AHE was initiated by Karplus and Luttinger back in 1954 
Q , and 20 year long efforts resulted in the conclusion that 
the AHE is an extrinsic effect controlled by competing 
terms related to electron scattering. Some of the contri- 
butions to Hall current found in this way do not depend 
of the scattering time 0, and recent research indicates 
existence of an intrinsic contribution to the anomalous 
Hall current that can be expressed in terms of Berry cur- 
vature in k -space 0,0,0- Experiment still does not pro- 



vide any convincing evidence regarding the mechanisms 
involved @, 0] . Identification of spin transport mecha- 
nisms is even a more challenging task because the very 
notion of spin current is not well defined in media with 
SO coupling. 



SPIN-ORBIT COUPLING 

In vacuum, SO coupling is described by the Thomas 
term in Pauli equation, 

H so { r ) = {eh 2 /Am 2 c 2 ) <x • (W( r ) x k ), (1) 

that is small for nonrelativistic momenta Tik -C tuqc , 
V( r ) being the scalar potential. In semiconductors with 
the gap Eg and SO band splitting A of a comparable 
magnitude, Eq ~ A ~ 1 eV, spin-orbit coupling is en- 
hanced by a factor of about moc 2 / Eq ■ The specific form 
of the SO contribution to the free electron Hamiltonian is 
controlled by the symmetry requirements 10]. Two typ- 
ical SO Hamiltonians for electrons in three-dimensional 
(3D) crystals are 

r(3D) 



H-r = a(er x k 



= P 3D (cr ■ K ), 



(2) 



where K z — k z (k 2 — k 2 ) , and K x and K y can be ob- 
tained from K z by cyclic permutations. The term H-r 
(Rashba term) is typical of hexagonal A 2 B g crystals. 
Being linear in k , it is the leading term of the k ■ p 
- theory. Hence, it makes a profound effect on electron 
spin dynamics [111. The term H^ 13 - 1 is typical of cu- 
bic A 3 B 5 crystals and is known as the Dresselhaus 3D 
term. 

In the Bloch function representation, SO coupling also 
changes the operator of coordinate r — iV ^ + r so . 
For electrons in narrow-gap semiconductors, r so equals 

r 90 = (/iV4mo)[£G 1 + (^G + Ar 1 ]((7 x k) (3) 

and is known as Yafet term g being the g -factor. 
In modern terms, r so is Berry curvature in k -space. 

Spatial confinement lowers the system symmetry and 
makes appearance of k -linear terms generic. In partic- 
ular, Hn develops because and is controlled by the con- 
finement asymmetry 0, or external strain 0] and 
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is known as the structure induced asymmetry (SI A) or 



Rashba term, while H, 



(3D) 



reduces to the bulk induced 



asymmetry (BIA) or Dresselhaus term 

Hb = /3(a x k x - <t y ky), (3 sw -/33D(7r/cf) 2 , 



(4) 



d being the confinement layer width. For heavy holes, 
k -linear terms are small and k 3 -terms dominate, 



hh 



/3hh0+c- - k 3 _a + ), 



(5) 



where k± — k x ± ife^ and er± = ct^ ± ia y |l6|. Relative 
magnitude of a and /3 is material and geometry depen- 
dent. It was measured by electrical 0,13' optical [lfij . 
and photoelectrical means. Interestin g eff ects of the 
interplay of -Hr and ifo were predicted |2lj, |22, |23| • 

The effect of SO coupling on spin transport critically 
depends on the analytical form of the SO Hamiltonian. 



INTERPLAY OF SPIN AND ORBIT: 

MAGNETIC SUSCEPTIBILITY 

Existence of SO interaction poses the questions (i) how 
deeply can this interaction influence macroscopic proper- 
ties of a coupled system, and (ii) under which conditions 
spin dynamics can be separated (at least approximately) 
from the orbital dynamics. It is natural to begin with an 
equilibrium property like magnetic succeptibility x ■ For 
2D systems, magnetization M strongly oscillates as a 
function of B , hence, the results are difficult to analyze 
[24|. Explicit results are available for the dependence of 
low- B susceptibility _v on the Fermi energy £p of 3D 
electrons at T = |25j. Hr electrons with isotropic 
mass m and g -factor in a magnetic field B || z are 
paramagnetic for small £f 



X 



e 2 (e F + s a ) 
Sitamc 2 



I - 



2ni{) 



e a < £f < 0; (6) 



here e a =ma 2 /h 2 . The susceptibility x vanishes when 
the Fermi level is at the spectrum bottom, £p = — £ Q , 
and increases linearly with ep ■ When the Fermi level 
passes through the singular point of the spectrum, £p = 
, the system becomes diamagnetic and x diverges 



X = 



in 2 he 2 \/2m£F ' 



< £ F < £ c 



(7) 



Further on, x( £ f) changes smoothly and asymptotically, 
for £f S> £ Q , approaches the Landau-Pauli limit. 

Therefore, for |ef| ~ e a spin and orbital motions are 
strongly coupled and their contributions cannot be sep- 
arated. It is only when £p 3* £a (and spin velocity 
a/h <C Dp ) that the orbital motion becomes a pertur- 
bation to spin dynamics, and long Dyakonov-Perel |2^| 
spin relaxation times can be achieved [!]]]• For interfer- 
ence SO devices exploiting exact eigenstates, fulfillment 
of this criterion is not required. 



BASIC IDEAS AND EXPERIMENTAL 
ACHIEVEMENTS 



A free-electron Hamiltonian with SO coupling Hr 
has, for each energy e , two eigenstates with the mo- 
menta k±(e) and opposite spin chirality. Datta and 
Das [3] initiated the idea of spin transistor based on 
injecting spin polarized electrons and the following in- 
terference of the two components of the incident beam. 
The conductivity of the device depends on the phase 
difference AO = (fc_ — k + )L acquired by these states 
after crossing a ballistic sample of a length L and os- 
cillates with a period defined by the interference condi- 
tion — k + )L = 2im , with n an integer. An equiv- 
alent description of the same phenomenon is the pre- 
cession of an electron spin in an effective magnetic field 
B so = (2a/gfj,s)( k xz) , [i^ being the Bohr magneton. 
The device is controlled by the gate voltage V s that mod- 
ulates the SO coupling constant a , a = ot(V g ) . This 
seminal paper enjoyed a lot of attention and generated 
active discussions and new proposals. I am not aware of 
any practical realization of this device. A diffusive analog 
of it with spin coherence restored with a spatial period L 
that does not depend on the driving electric field E has 
been demonstrated recently 29] ; this periodicity proves 
that SO coupling was linear in k . 

Koga et al. [30] proposed a ballistic spin-interfcrometer 
using square loop geometry with the interference of par- 
tial waves at the incident point. It does not require inject- 
ing spin-polarized electrons, is free of magnetic elements, 
and is completely controlled by the gate voltage V g . Re- 
markably, they achieved a phase difference AO ~ tt per 
loop for an extensive loop array; conductivity modula- 
tion was detected by the V„ dependence of Al'tshuler- 
Aronov-Spivak oscillations |31| . 

In an external magnetic field, spin polarized ballistic 
electron beams were created due to SO coupling in an 
open quantum dot serving as an emitter |32| or in the 
bulk [33j. For a spin-split energy spectrum, refraction 
and reflection of a spin-unpolarizcd electron beam trans- 
forms it into a system of spin-polarized beams |34l l3"J| , 
and spin polarization of electrons scattered from a litho- 
graphic barrier has been achieved [36|. 

Recently, experimental discovery of spin-Hall effect 
(SHE) with 3D electrons H3 and 2D holes HU attracted 
attention. Dyakonov and Perel |3^] predicted SHE due to 
Mott' skew scattering in a diffusive regime. It manifests 
itself in spin accumulation near the sample edges even 
in the absence of an external magnetic field. Because 
shortly before this discovery dissipationless spin currents 
for 3D holes in diamond type materials ^(J an d f° r 2D 
electrons in systems with SI A have been proposed, 
the mechanism of SHE became a matter of active dis- 
cussions. Below, I discuss this problem as applied to 2D 
systems. 
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SPIN CURRENTS AND MAGNETIZATION 

A physically measurable quantity is magnetization, 
and when SO coupling is not too strong, £p ^> e a , one 
can distinguish orbital and spin contributions to it. It is 
appealing to discuss spin accumulation in terms of spin 
currents generated due to SO interaction and driven by 
an electric field E ; such an approach was employed by 
Governale et al. [42| and Mal'shukov et al. |43J. In 
vacuum, SO coupling exists only near scatterers and is 
described by Eq. Scattered electrons are polarized 
along ( k x k ') , where fc and k ' are the initial and 
final momenta, respectively. Electrons scattered to the 
left and right have opposite polarizations. Spin currents 
are conserved in the free space. 

In semiconductors, in addition to this extrinsic effect, 
there exists also an intrinsic SO effect originating due 
to the spin dependence of the free electron Hamiltonian 
H(k ) . Calculating the cumulative effect of both mecha- 
nisms to spin transport in terms of spin currents is prob- 
lematic for several reasons. 

(i) Because of the torque produced by the field B so , 
spin is not conserved and does not obey the continuity 
equation. Therefore, the standard procedure cannot be 
applied for defining spin currents. Usually they are de- 
fined similarly to the media without SO coupling as 

J li = \{ v i a i +Vj v i)i i = {x,y}, j = {x,y,z}, (8) 

where the velocity v is defined, with the due account of 
the SO coupling in the bulk, as v = h dH{ k) Id k . 
Competing definitions of J ] si were also proposed j44j. 

(ii) Spin currents of Eq. © are even with respect to 
time inversion, t — *■ — t , hence, some of the components 
of J 3 si may be nonzero in equilibrium. E.g., for the 
Hamiltonian Hr , J| x = — J* y ^ 0. In this respect, 
spin currents are reminiscent of the gas pressure that is 
the momentum flux. Similarly to it, the currents J 3 si , 
depending on the conditions, can play a role of both ther- 
modynamic and transport parameters |45| . 

(iii) Because in equilibrium J° s i ^ in media with SO 

coupling, and J 3 S i — in media without it, spin currents 
are not conserved at their interfaces. 

(iv) J 3 si are i-even while the magnetization M is 
t -odd, hence, magnetization build-up near a boundary 
should involve some t -inversion breaking process. 

Because of these problems, bulk spin currents of Eq. (JSJ 
seem not be directly related to spin accumulation near 
edges. Nevertheless, calculations of them turned to be 
highly instructive regarding the dependence of spin trans- 
port on SO coupling mechanisms. 



EFFECT OF ELECTRON SCATTERING ON 
SPIN CURRENTS 

Dissipationless spin currents were found originally by 
applying the Kubo formula in the regime when scattering 
was disregarded. Such an approach is questionable be- 
cause the linear response provided by this formula can 
be maintained only due to scattering. A proper ac- 
count for electron scattering resulted in different results 
for the Hamiltonians and , and the Luttinger 
Hamiltonian of 3D holes. For Hn electrons, Inoue et al. 
[lif , followed a number of different authors 0, EM Oil > 
have found dramatic cancelation of the bubble and lad- 
der diagrams for J* x , hence, spin current vanishes, and 
this is true for any type of elastic scattering |50|. On 
the contrary, for more general SO coupling mechanisms, 
Jsx 7^ |5lL 152. l53| . In particular, for the Hamiltonian 
ifhh ladder diagrams vanish [HjJ, and J* x is described 
by bubble diagrams 0- 

Dimitrova 0] related spin-current cancelation for the 
operator to the operator identities 

fr x = -k a J* x , fry = -k a Jl y , k a = am/h 2 , (9) 

k a being the spin-precession wave vector. Eq. @ im- 
plies that non-zero macroscopic spin currents, J s z x =/= , 
Jay > would be tantamount to a time dependent 
macroscopic magnetization. This argument is valid for 
a generic combination of and Hd (a ^= ±/3); for 
more detail see Refs. [5^ and [5(|. It also persists for 
spin-independent electron scattering. 

A different insight into spin-current cancelation comes 
from applying the Kubo formula to electrons in a 
mag netic field B || z (5^|. This system is exactly soluble 
, and for calculating transverse responses introducing 
impurity scattering is not required. Spin current J* x 
includes two contributions. First comes from interbranch 
transitions and includes a sum over the Fermi sea region 
with fc + < k < k- . Second comes from two intrabranch 
transitions at the Fermi sea edges, one per branch. For 
g = , these terms cancel for arbitrary B and a/hvF 
because of exact sum rules specific for Hamiltonian, 
hence, J* x — . They ensure insensitivity of J s z x to 
scattering mechanisms. 

For B — > , the separation Tilo c between adja- 
cent levels vanishes, u c being the cyclotron frequency. 
Thence, the intrabranch contribution to spin conductiv- 
ity is highly sensitive to any inaccuracies in numerical 
simulations, while the interbranch contribution equal to 
the universal conductivity ct s h = e/Airh by Sinova et al. 
411 is robust. That is why special efforts were needed 
53. l6Cl| to bring numerical results on Hamiltonian 
in agreement with the analytical ones. 
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SHE: FREQUENCY AND MOMENTUM 
DEPENDENCE 



Sensitivity of intrabranch terms to the external param- 
eters manifests itself in a strong frequency dependence at 
the scale of w ~ uj c 58] . With oj — > i/r , where r is the 
momentum relaxation time, ct s h(w) matches the result 
of Ref. 0]- I n the frequency range uj c <C ui <C 2ak±/h , 
spin conductivity <t s h(w) « e/Airh. 

Similarly, er S H shows momentum dependence at the 
scale of k ~ A: a |45j. To find spatial dependence of 
spin accumulation S ( r ) , we take advantage of using 
diffusive equations derived by Mishchenko et al. 0] an d 
Burkov et al. For i/^ electrons in the x > half- 
plane driven by a field E || y , 

( j D^'-5 x /t s ) + 2Mvf^ = S£7t s , 
(DS" - - 2k a £v F S' x = 0, S y = 0, (10) 



where £ = vpr , 13 



! p t/2, r s = l/2r(fc Q w f ) 2 , and 



S!j? = k a TeE/2ir are the mean free path, diffusion co- 
efficient, spin relaxation time, and spin accumulation at 
x — > oo , respectively. has been measured recently 

Eq. (HOI is valid when k a l < 1 . The boundary 
conditions for it should relate S and S ' at x = . 



SHE due to spin-current cancelation breaking should 
also develop in systems with non-planar boundaries and 
inhomogeneities. Stern-Gerlach s pin filter due to a = 
a(x) was proposed by Ohe et al. |57j . 

Remarkably, SHE of Fig. 1 has been found for iJ K 
electrons, i.e., for a system with J* x = . This fact in- 
dicates that bulk currents Jf ■ are not enough for calcu- 
lating SHE, and transport equations valid at the spatial 
scale of £ a — k~ 1 are needed. For iJhh holes they are 
not available yet, hence , consistent treatment of the data 
by Wunderlich et al. is not possible. However, if one 
accepts that transport spin currents at the spin preces- 
sion wave vector fc so = kp ~ m^hhk^/h are governed 
by "macroscopic" spin conductivity and are about eE /h 
[53| . one can speculate about SHE. Indeed, S z (0) ~ 
hJ^ x /v e g , and v e s ~ k^ 1 Jt is a proper estimate for 
spin-transport velocity. Hence, ^(O) ~ fc^reS, simi- 
larly to found above, and large SHE of Ref. [38| can 
be attributed to large magnitude of kpr . 

Generalizing the above conclusions, we make a Con- 
jecture that, irrespective to specific SO coupling mecha- 
nisms, spin currents at the wave vector k so have mag- 
nitudes about the "universal" value of eE/Ti 0] and 
acquire the meaning of transport currents. Hence 



JsH(fcso) ~ eE/Ti , S(k so )/h ~ k so reE/h . 



(11) 
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FIG. 1: Spin accumulation of Hb. electrons in the x > 
half-plane driven by a field E || y . Upper curves: S x , 
lower curves: S z ; all in units of S£° . Dimensionless spin 
relaxation velocity s = s/k a D , the coordinate x in units of 



In Fig. 1 are shown the results for the boundary con- 
dition Z)S-(O) = sSi(0) , i = x,y, where s is the 
surface spin-relaxation velocity. The condition s « 
sounds reasonable for smooth lithographic barriers along 

free edges, hence, we corroborate the conclusion by 

Mishchenko et al. [47j regarding the absence of spin accu- 
mulation along them. By engineering high s edges, one 
can achieve spin accumulation along them. All curves 
in Fig. 1 are scaled in units of , hence, the process 
looks like rotation of S from the x to z direction due 
to the transverse spin flow. Leaking of spins into a = 
reservoirs results in a similar effect 56] . 



Numerical coefficients in these equations depend on the 
specific form of H so and on boundary conditions. 

Difference in SHE in meso- and macroscopic Hr con- 
ductors, emphasized in Ref. jrJij ]. seems to come from the 
momentum dependence of spin conductivity. It also en- 
ables spin injection from quantum rings |65, 66]. 



EXTRINSIC SHE IN JV-GAAS 

Engel, Halperin and the present author [67j have de- 
veloped a theory of extrinsic SHE in n -GaAs. Its mag- 
nitude comes from the enhancement of the coefficient in 
Eq. JQ by a factor of more than 10 6 due to the cou- 
pling across the gap Eq . A Boltzmann theory of SHE 
results in skew scattering and side jump terms that are of 
comparable magnitude for materials with small k-pt m 4 
used in Ref. |37| . Remarkably, the side jump term in- 
cludes a Berry phase part that might indicate a way for 
unifying the intrinsic and extrinsic approaches to SHE. 
The results are in a reasonable agreement with the data 
by Kato et al. piTl. More recently, Tse and Das Sarma 
came to similar conclusions 681. 
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